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Abstract. In order to investigate the nature and origin 
of dust in elliptical galaxies, a method for the dust mass 
evaluation, which accounts for the dust temperature dis- 
tribution, is here presented and discussed. The derived 
dust masses turn out to be a factor 2-6 larger than those 
obtained with the single temperature approximation. A 
correlation between the far-infrared and the blue luminos- 
ity has been also found. The results are discussed in terms 
of dust "mass discrepancy" and of the possible evolution 
scenarios: evaporation flow and/or cooling flow. While the 
present data cannot discriminate between these two sce- 
narios, it is conceivable that the dust in elliptical galaxies 
can be accreted by the contribution of different mecha- 
nisms, according to the history of the individual objects. 

Key words: Galaxies: general - Galaxies: elliptical and 
lenticular, cD - Galaxies: ISM - Cooling flows - Infrared: 
galaxies - Infrared: ISM: continuum 



1. Introduction 

Evidence for the presence of dust in elliptical galaxies was 
given by the optical observations of obscured regions in 
some systems (see for instance Bertola et al. 1985, Veron 
& Veron 1988) and was finally confirmed by the FIR obser- 
vations of the IRAS satellite (Knapp et al. 1989, Roberts 
et al. 1991). The low resolution IRAS data (> 1 arcmin) 
have to be processed by adequate techniques to provide 
detailed (~ 1 arcmin at 100 /im) information about the 
dust spatial distribution in nearby elliptical galaxies. In 
general, only the integrated flux of the detected source is 
available in the IRAS bands. 

Despite the intrinsic limits due to the low resolution, 
Knapp et al. (1989) showed that a significant fraction 
(48%) of the nearby E and SO galaxies from the Revised 
Shapley-Ames Catalogue (Sandage & Tammann 1981, here- 
after RSA) have been detected by IRAS at 60 and 100 /im 



at the limiting sensitivity (about 3 times lower than in the 
IRAS Point Source Catalog). 

It is not surprising that ellipticals contain dust, since 
the presence of dust is directly related to the stellar forma- 
tion. But the coexistence of solid particles with the dom- 
inant gas component, which in these galaxies is heated to 
~ 10 7 K and radiates at X-ray wavelengths, is a matter 
of discussion. In fact, dust grains should be quickly de- 
stroyed by sputtering (Draine & Salpeter 1979) when in 
direct contact with the hot gas. In such an environment 
the dust has a lifetime of 10 6 — 10 7 yr. The question thus 
follows: where does the dust come from? 

At FIR wavelengths (60 and 100 /im) the thermal emis- 
sion is mainly due to large grains (radius ~ 0.1 /im) which 
may have different heating sources, e.g. the general inter- 
stellar radiation field or OB stars. In order to discriminate 
between the two different contributions and to estimate 
the weight of each of them, several efforts were under- 
taken (see for instance Calzetti et al. 1995). A reasonable 
approach is to consider the 60 /im flux entirely due to the 
warm dust (~40 K), while the 100 /im flux should be con- 
sidered the result of two contributions: the warm and the 
cold (^10 K) dust. While a two-component dust model 
is used to explain the spectral trend for different types of 
galaxies, it is rarely adopted for elliptical galaxies, which 
are often characterized by weak FIR emission and which 
are not always detected in all IRAS bands. Therefore, a 
single color temperature (from the 60 and 100 /im data) 
is usually taken as the dust temperature. It follows that 
no information about the dust temperature distribution 
and the dust spatial distribution is available for elliptical 
galaxies. The availability of the ISO data will provide spec- 
tra in a wider IR wavelength range (2.5-240 /im). Several 
attempts to understand the dust nature and origin suggest 
interesting interpretations by comparing optical and FIR 
data (Goudfrooij & de Jong 1995, hereafter GJ95, Tsai & 
Mathews 1995, 1996), by studying the stellar content, or 
by using a severe and critical approach to the data (Breg- 
man et al. 1998). Finally, few elliptical galaxies have been 
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observed at sub-millimeter wavelengths by Fich & Hodge 
(1991, 1993). 

GJ95 found that the dust masses determined from the 
IRAS flux densities are roughly an order of magnitude 
higher than those determined from optical extinction val- 
ues of dust lanes and patches, in contrast with what hap- 
pens for spiral galaxies. The authors suggest that this 
"mass discrepancy" may be explained by the existence 
of a diffuse component (within 2 Kpc from the center), 
which is not detectable at optical wavelengths. On the 
other hand, Tsai & Mathews (1996) suggested that, while 
the distributed dust component is associated with dust re- 
cently ejected from evolving stars, another "extra" compo- 
nent of dust is present in ellipticals both in dust lanes and 
rings and/or in other galactic regions. In particular, they 
postulated that a substantial mass of cold gas remains 
"observationally elusive without forming completely into 
stars" . If the extra dust is optically thin in the visible it 
should be located far from the galactic core region, where 
the intensity of the starlight and therefore the grain heat- 
ing is reduced. 

The dust spatial distribution suggested by GJ95 and 
by Tsai & Mathews (1996) support the two most popu- 
lar scenarios respectively: the evaporation flow picture and 
the cooling flow picture. In the evaporation flow scenario 
the clouds of dust and gas currently observed in ellipticals 
have mainly an external origin, being associated to events 
of galaxy interaction and/or mergers and being heated by 
thermal conduction in the hot gas (de Jong et al. 1990, 
Sparks et al. 1989). On the contrary (cooling flow), the in- 
ternal origin of gas and dust may be explained with both 
red giant winds (Knapp et al. 1992) and by the cooling 
flow mechanism, in which mass loss from stars within the 
galaxy is heated by supernova explosions and by collisions 
between expanding stellar envelopes during the galaxy for- 
mation stage, and then cools and condenses (Fabian et 
al.1991). 

Therefore, in ellipticals, the presence of dust and the 
dust "mass discrepancy" are related with the dust spatial 
distribution, which depends on the nature and the evolu- 
tion of these systems. Hence, in order to investigate the 
dust content, dust mass evaluations as accurate as pos- 
sible are required to estimate the amount of the "mass 
discrepancy". Since both scenarios suggest the presence 
of different dust components and the very nature of the 
galactic environments involves the existence of dust grains 
at different temperatures, it is necessary to take into ac- 
count a dust temperature distribution for the dust mass 
evaluations. 

Kwan & Xie (1992) suggest a theoretical approach to 
take into account the effects of the dust temperature dis- 
tribution in the dust mass evaluation. I present here an 
application of their method which is discussed and imple- 
mented in Section 2. The results obtained for the sample 
of ellipticals, introduced in Section 3, are presented and 
discussed in Section 4 by comparing the FIR with the 



visual dust mass evaluation and by discussing the corre- 
lation between blue and FIR luminosities. 



2. Dust mass evaluation from FIR observations 

The dust mass can be derived both from optical and from 
FIR observations. The value of the mass depends on the 
physical-chemical properties of the solid particles (i.e. grain 
radius, grain density and emissivity). The usual approach 
consists of assuming some values for the grain proper- 
ties (Hildebrand 1983), while a color temperature is de- 
rived from the FIR emission. Different authors, however, 
present slightly different formulae for the evaluation of the 
dust mass (Thronson & Telesco 1986, Greenhouse et al. 
1988, Young et al. 1989, Roberts et al. 1991 and Thuan & 
Sauvage 1992). The differences between these relations are 
only due to different assumptions on the grain parameters 
and to different derivations of the color temperature. It 
should be noticed that, within a flux uncertainty of 10% 
(which is typical of the high quality IRAS data), all the 
dust mass values obtained by using the different formulae 
are in agreement. 

A single temperature model is a rough approximation 
in describing a galactic environment. The dust is in fact 
heated by the radiation field, which in turn depends on the 
sources of luminosity and on their spatial distribution in 
the galaxy. The total FIR emission is thus likely due to the 
contribution of dust at different temperatures. Moreover, 
the IRAS FIR measurements are not adequate to detect 
the emission coming from cold dust (10 — 20 K) which 
peaks at wavelengths between 200 and 300 /im. 

I adopt the dust temperature distribution given by 
Kwan & Xie (1992): 

(T/TrVe-^ T / TL ^ 
9(T) = , T }, al „,„ ; f° r Tl <T < Tjj , (1) 



J^{T/T L pe-f 3 ( T / T ^)dT 

g(T) = otherwise. (2) 

Tl and Tjj are the lower and upper limits of the temper- 
ature T; [3 and 7 are free parameters that determine the 
shape of the distribution. The equations relating the tem- 
perature distribution g(T) to the luminosity emitted by 
the dust and to the dust mass are detailed in Kwan & Xie 
(1992). 

Due to the observed spectral range, I adopt Tj,=l K 
and T;y=60 K, taking into account only temperature dis- 
tributions peaking at a value intermediate between Tj, and 
Tjj and excluding those pairs of j3 and 7 which produce 
unrealistic distributions (as, for instance, monotonically 
increasing or decreasing functions). 

The main problem in the present approach is to select 
the proper values for the parameters (3 and 7, the choice 
being constrained by the ratio of the flux densities at two 
different wavelengths. 

I first identify a range of [3 and 7 pairs which pro- 
duce the observed flux ratio. Since the same flux ratio can 
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be obtained by functions having quite different shapes, a 
further constraint is needed. Unfortunately, for none of 
the galaxies in the sample submillimcter observations are 
available. Nevertheless, their color temperature may be 
derived from the flux density ratio (Hcnning et al. 1990). 
Taking into account the uncertainty in the computed color 
temperature, I select those distributions g(T) that pro- 
duce the observed flux density ratio and whose peak tem- 
perature is comparable to the color temperature. The se- 
lected family of g(T) functions obviously satisfies two con- 
ditions which arc not independent (being both related to 
the flux ratio). This fact could in principle affect the reli- 
ability of the results. It turns out, however, that pairs of 
(3 and 7 which satisfy the same constraints give the same 
dust mass within the flux uncertainty. 

In order to check the method a simple numerical simu- 
lation has been performed. I considered an artificial galaxy 
which is not resolved by IRAS both at 60 and 100 ^m with 
a given dust mass and, then, I evaluated the dust content 
following the present method. The galaxy dust mass Md 
may be roughly estimated by using the equation 

M d = ^r 3 P N d , (3) 

where r is the dust grain radius, p is the dust grain density 
and Nd is the number of dust grains. Eq. 3 implies the fol- 
lowing approximations. The dust grain radius r character- 
izes the dust grains which contribute to the thermal emis- 
sion at FIR wavelengths. Actually, the uncertainties in the 
total dust mass introduced by using an average radius in- 
stead of a size distribution, are much smaller than those 
arising from uncertainties in the dust emissivity spectral 
trend and in the dust temperature distribution (Kwan & 
Xie 1992). is the total number of the dust grains within 
the observing beam. Since the galaxy is not resolved by 
IRAS, as it is often the case, Eq. 3 gives the total dust 
mass of the point source. A dust radius of 0.1 /im and 
a dust grain density of 3 g cm -3 are currently adopted 
(Hildebrand 1983) and also used in the present article. 

Concerning the dust emissivity, the power-law approxi- 
mation e(A) = eo(Ao/A) a is used in the computations, with 
a = 1 and e = 9.38 • 10~ 4 (Hildebrand 1983). 

By assuming a value for N d and the temperature dis- 
tribution g(T), and accounting for the spectral responses 
of the IRAS detectors at 60 and 100 /im, the luminosity 
emitted by the Nd dust grains, as observed within the fil- 
ter bandpass, can be estimated. The flux ratio and the 
color temperature of the galaxy are then computed, while 
the dust mass of the source is derived from Nd by Eq. 3. 

A galaxy with Md=10 5 M Q turns out to contain about 
A,i=1.6-10 52 dust grains. The galaxy distance is assumed 
to be 10 Mpc. By assuming a g(T) with (3 = 5 and 7 = 20 
the flux ratio turns out to be about 0.5. I use the flux 
ratio as a constraint to select the pairs /3-7 producing a 
g(T) whose peak temperature is comparable to the galaxy 
color temperature of 36 K. Among these pairs also the pair 
P = 5 and 7 = 20 is found (i.e. exactly the values adopted 



Table 1. The Selected Sample of Elliptical Galaxies. 
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Column (2): total blue luminosities taken from GJ95. Columns 
(3), (4) and (5): computed values for FIR luminosities, flux 
ratio and dust color temperatures. The dust masses in Column 
(6) are taken from GJ95. The FIR dust mass in Column (7) is 
derived by averaging the results of formulae which are usually 
adopted, while the FIR dust mass in Column (8) is derived 
by adopting a temperature distribution model as described in 
Sect. 2. The uncertainties for the mass evaluations are shown 
in Fig. 1. N and I stand for NGC and IC respectively. 

for the input g(T)). For the selected pairs of parameters I 
compute the dust mass which ranges between 6 — 9 • 10 4 
Mq. Therefore, the derived values are in agreement with 
the dust mass of the galaxy within the flux uncertainties, 
thus confirming the reliability of the method. Futhermore, 
by using the single temperature model and taking into ac- 
count the different formulae available (see sect. 2) a dust 
mass of 4 • 10 4 M Q is derived. This result shows that the 
single temperature model underestimates the dust con- 
tent. 

The same test was performed for different "artificial 
galaxies" with different temperature distributions obtain- 
ing always consistent results. 

3. The sample of elliptical galaxies 

The 21 elliptical galaxies listed in Table 1 have been ex- 
tracted from the sample by GJ95, who evaluated their FIR 
luminosity and dust mass from FIR data. All the galaxies 
are classified as E both in RSA and in de Vaucouleurs et 
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al. (1991) and have magnitudes < 12. I selected those 
galaxies which have been detected at both 60 and 100 /xm 
by IRAS (Knapp et al. 1989). H = 50 Km s^Mpc" 1 is 
assumed throughout this paper. 

For each galaxy the FIR luminosity the flux ratio, the 
dust temperature and the dust mass from FIR data (here- 
after FIR dust mass) have been derived and listed in Ta- 
ble 1. The FIR fluxes by Knapp et al. (1989) have been 
corrected taking into account the contribution of hot cir- 
cumstellar dust (GJ95) and have been used for the com- 
putation of: the FIR luminosity in the band 40-120 /im 
(Helou et al. 1985), the flux ratio and the dust color tem- 
perature (Hcnning et al. 1990). The uncertainties in the 
three quantities depend on the flux uncertainties, which 
are estimated to be in the range 10%-30% (Knapp et al. 
1989). The chemical and physical properties of the dust 
grains play a critical role in the evaluation of the uncer- 
tainties on temperature and mass. On the other hand, 
the present paper is mainly addressed to the dependence 
of the dust mass on the choice of a particular tempera- 
ture distribution, and I will not discuss the problem of 
the dust properties (which will be studied in a forthcom- 
ing paper). A spectral index a = 1 for the dust is used in 
all the computations. The FIR dust masses listed in Table 
1 are evaluated from the FIR fluxes by means of the sin- 
gle temperature and the temperature distribution models 
respectively. 
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Fig. 1. The dust masses computed by adopting a single 
temperature (open squares) and a temperature distribu- 
tion (filled circles). The error bars account for the uncer- 
tainties in g(T) (filled circles) and in the color temperature 
(open squares). 



4. Temperature distribution model: results and 
discussion 

4-1. FIR dust masses 

The dust masses computed by adopting a single tempera- 
ture and a temperature distribution are compared in Fig.l. 
Both methods lead to a temperature dependence, as ex- 
pected when the dust amount is derived by using thermal 
emission. In fact, the colder the dust grains, the larger 
should be their total number in order to produce a given 
FIR emission. 

The dust masses obtained with a temperature distribu- 
tion are larger than those derived from a single-tempera- 
ture model. This confirms that IRAS measurements allow 
to determine the warm (30-54 K) dust amount only, while 
the contribution of the cold dust is neglected. Moreover, it 
turns out that, in the present sample, the FIR luminosity 
does not depend on the color temperature, thus confirming 
that the FIR emission is a result of different contributions 
and cannot be properly explained by a source with a single 
equilibrium temperature. The two models give masses dif- 
fering by factors from 2 to 6. This wide range is due either 
to the shape of the temperature distribution and/or to the 
uncertainties in the dust parameters. On the other hand, 
the ratio between the two dust mass evaluations (Fig. 2) 
shows a general temperature dependence which suggests 
that colder galaxies have a larger amount of missed dust. 




30 40 

Dust Color Temperature [K] 



50 



Fig. 2. The ratio between the FIR dust masses computed 
by using a temperature distribution and a single temper- 
ature model. The temperature dependence suggests that 
larger amounts of cold dust are neglected in sources with 
lower color temperature. The two galaxies on the left bot- 
tom, which do not follow the general trend are, from left 
to right, NGC 1395 and NGC4589. 
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4-2. Dust "mass discrepancy" 

The temperature distribution model enhances the dust 
"mass discrepancy" . In order to explain that discrepancy, 
it has to be noticed that the dust mass evaluated from 
optical observations critically depends on the spatial dis- 
tribution of the dust with respect to the stars. Usually the 
extinction is thought to be due to an overlying absorbing 
screen of dust grains. That is the geometry where a given 
amount of dust has the strongest effect on the starlight. 
In order to test how a peculiar and oversimplified spatial 
distribution affects the dust mass evaluation from opti- 
cal data, I compared the optical dust masses in Table 1 to 
those derived by assuming a more realistic spatial distribu- 
tion (e.g. Witt et al. 1992) and by using the optical data by 
Goudfrooij et al. (1994a, b). The new optical masses turn 
out to be enhanced by a factor 2-4, cutting the "mass dis- 
crepancy" down. However, even if the "mass discrepancy" 
is reduced by assuming a more realistic spatial distribu- 
tion, the optical absorption can only account for the dust 
in the obscured galaxy regions, while the distributed dust 
component is always neglected. 

A similar case concerns the "extra" dust component 
suggested by Tsai & Mathews (1996) to explain the 60- 
100 /Ltm flux ratio: due to its low temperature, it is not 
detected by IRAS. Following their model, it is possible 
to estimate the contribution of the "extra" dust in the 
sample. Taking into account the "extra" dust amount, the 
dust masses derived by a single temperature model are 
slightly enhanced, but the large uncertainties due to the 
observations and, then, to the temperature evaluation do 
not allow to judge this enhancement as significative. One 
can argue that, with a model of temperature distribution, 
it is possible to overestimate the dust amount. Actually, 
this cannot happen because of the severe constraints which 
are chosen, whether two different flux ratios are available 
or when a flux ratio and the relative color temperature 
are used as it is here suggested. In this respect, the last 
method is also the more conservative, since the warm dust 
observed at 60 /im affects the dust mass evaluation in the 
sense of underestimating the cold dust amount. 

Therefore, since these computations are very sensitive 
to the assumed dust temperature, the model by Kwan & 
Xic (1992) greatly improves the mass evaluation by taking 
into account a wide temperature range, and can be applied 
to a wide range of galaxy morphological types when, be- 
cause of lack of the data, it is no possible to adopt and to 
develop a suitable radiation model. 

5. Temperature distribution and heating 
mechanisms: the origin of the dust 

Actually, the uncertainty on the temperature distribution 
is directly connected to the study of the heating mecha- 
nisms and, therefore, to the knowledge of the UV inter- 
stellar radiation field. 
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Fig. 3. a) Dust mass (with the temperature distribu- 
tion model) and blue luminosity relation: the solid line 
is the least squares bisector fit (confidence level greater 
than 95%) while the dashed line is the loci where dust is 
replenished by stellar mass loss and destroyed by sputter- 
ing with the maximum destruction timescale (10 yr). 
Most of the galaxies present an observed dust mass larger 
than the amount predicted by the dust formation and de- 
struction model, b) FIR and blue luminosity relation: the 
solid line is the least squares bisector fit (confidence level 
greater than 99.5%). For comparison the theoretical pre- 
diction (dashed line) by Tsai & Mathews (1996) is also 
shown. 

To this aim, and in order to investigate the origin of 
the dust in elliptical galaxies, I plot in Fig. 3 a and b the 
dust mass M<j and the FIR luminosity Ijfir versus the 
blue luminosity Lb for the galaxy sample. Both M<j and 
L fir are correlated to Lb , in contrast with what has been 
found by GJ95. Therefore, while the absence of correlation 
between dust mass and blue luminosity in their sample 
was considered to support the external origin of the dust 
in ellipticals and the evaporation flow picture, the present 
result indicates that a relationship between the dust con- 
tent and the present day population of stars should be not 
excluded. But, as already pointed out by GJ95, a lack or a 
presence of correlation between M^ and/or Lfir and Lb 
cannot be a definitive proof of the external or the internal 
origin of the dust, due to the fact that the dust destruc- 
tion mechanisms and timescales in elliptical galaxies are 
expected to be different depending on the evolutionary 
state and on the hot gas content of the individual objects. 

For this reason, further analysis and, above all, high 
resolution FIR observations are required to understand 
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the origin and the fate of dust in these systems. Waiting for 
a complete analysis of the ISO data, a reasonable approach 
is the comparison between the available observations and 
models of star formation and dust heating mechanisms. I 
plot in Fig. 3a the loci where dust is replenished by stellar 
mass loss at the rate given by Faber & Gallagher (1976) 
and destroyed by sputtering with the maximum destruc- 
tion timescale (10 7 5 yr) in ellipticals which contain hot 
gas (see GJ95 for details). In 85% of the galaxies the dust 
having an internal origin does not account for the com- 
puted dust masses, which are larger. Therefore, one can 
argue that in 85% of the ellipticals an alternative supply- 
ing mechanism for the dust is required to account for the 
FIR observations. Mergers with spirals or small, dust-rich 
irregular galaxies and , then, the evaporation flow scenario 
are, in fact, strongly supported by the previous evidence 
(Fig. 3a). The critical point of this picture is to identify 
the observed diffuse dust, located within 2 Kpc from the 
center, as the dust accreted during the galaxy interac- 
tion when high hydrodynamic instabilities are expected. 
In particular, the external dust had to be somehow pro- 
tected from the interaction with the hot gas when moving 
toward the center. 



The problem could be solved with the "extra" dust 
component (Tsai & Mathews 1996) located in very large 
disks out to the effective radius and, therefore, cold enough 
to emit at FIR and submillimeter wavelengths. This dust 
can have, in principle, both internal and external origin. In 
particular, due to the low temperature "dust may re-form 
and grow in these cold disks" (Tsai & Mathews 1996). Un- 
fortunately, due to the large uncertainties, it is not possi- 
ble to evaluate the exact amount of this dust component. 
Therefore, since the external origin of the dust is far from 
fully verified, I experimented with an alternative approach 
to interpret the relation between the FIR and the blue lu- 
minosity. By using the least squares bisector method, I 
find Jjfir ocL^ 64±0 ' 39 (Fig. 3b), in agreement with Breg- 
man et al. (1998). The luminosity correlation supports the 
scenarios in which a significative amount of dust has in- 
ternal origin (Tsai & Mathews 1995, 1996), coming from 
stellar mass loss and being heated by stellar photons and 
the general interstellar radiation field. 

The problem of the coexistence of dust and hot gas 
might be resolved with a dust grain distribution contain- 
ing grains larger than the maximum size suggested by 
Mathis et al. (1977) and usually adopted in the current 
models. In fact, this assumption would increase the sput- 
tering time and, then, the dust grain density. Tsai & Math- 
ews (1996) found that the FIR luminosity is proportional 
to hg 1 when the maximum grain size increase from 0.3 
/im to 0.9 /im. This theoretical prediction (dashed line in 
Fig. 3b) and the observed relation (solid line in Fig. 3b) is 
consistent in the luminosity range covered by the present 



sample^, therefore the luminosity correlation seems to sup- 
port the cooling flow scenario. 

Finally, the two comparisons between observations and 
empirical-theoretical models show a clear conflict: while 
the canonical star formation rate cannot account for the 
dust amount in elliptical galaxies, thus supporting the 
evaporation flow scenario, the trend of the hpiR-^B can 
be explained with a proper dust model which increases 
the sputtering time and which suggests an ad hoc dust 
distribution. Concerning this last possibility it has to be 
stressed that, although the L^y^-Lb relation can be af- 
fected by different hot gas contents and stellar popula- 
tions in the individual objects, the correlation between 
these quantities is characterized by a trend predicted by 
theoretical models. Furthermore, both in Bregman et al. 
(1998) and in the present work, severe selection criteria 
are applied (see Sect. 3), which take into account the reli- 
ability of the observations and the different contributions 
to the FIR emission. I compare the present sample and 
the sample of Bregman et al. (1998) which rejected the 
galaxies whose FIR fluxes can be contaminated by AGN 
emission, position uncertainties, background objects and 
inhomogeneity. The Jjfir-Lb relation in Fig. 3b is not 
affected by the exclusion from the sample of the four 
AGNs (NGC4374, NGC4486, 14296 and 11459) and of 
NGC 3557 that has a non-homogenous background. 

Taking into account the results of the present analysis 
it seems difficult to choose between the evaporation flow 
or the cooling flow scenario and to affirm that the dust 
in ellipticals has mainly en external or an internal origin. 
A resonable conclusion is that the individual objects have 
experienced different mechanisms of dust accretion dur- 
ing their evolution and that only a detailed study of the 
individual systems will allow to describe for each object 
the different evolution phases and the contributions of the 
different mechanisms of dust accretion. 
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